Time-optimal Unitary Operations in Ising Chains II: 
Unequal Couplings and Fixed Fidelity 
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We analytically determine the minimal time and the optimal control laws required for the real- 
ization, up to an assigned fidelity and with a fixed energy available, of entangling quantum gates 
(CNOT) between indirectly coupled qubits of a trilinear Ising chain. The control is coherent and 
open loop, and it is represented by a local and continuous magnetic field acting on the intermediate 
qubit. The time cost of this local quantum operation is not restricted to be zero. When the match- 
ing with the target gate is perfect (fidelity equal to one) we provide exact solutions for the case of 
equal Ising coupling. For the more general case when some error is tolerated (fidelity smaller than 
one) we give perturbative solutions for unequal couplings. Comparison with previous numerical 
solutions for the minimal time to generate the same gates with the same Ising Hamiltonian but with 
instantaneous local controls shows that the latter are not time-optimal. 

PACS numbers: 03.67.-a, 03.67.Lx, 03.65.Ca, 02.30.Xx, 02.30.Yy 
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I. INTRODUCTION 

Time-optimal control is an active branch of the wider 
arena of quantum control. Quantum control has been 
applied successfully to a series of physical problems, 
e^., molecular dynamics [l[, quantum information 0]- 
[3(, multidimensional advanced spectroscopy techniques 
in nuclear magnetic resonance (NMR) experiments 
Time-optimal quantum computing aims at minimizing 
the time cost to achieve a certain target. It may be of 
relevance both from the theoretical point of view, e.g., 
in order to determine the fundamental maximal speed at 
which quantum information may be transported [H-Q, 
to find out faster quantum algorithms and solve com- 
plex problems Q, or it may serve more practical but 
essential purposes such as determining the fastest way to 
realize a certain unitary gate before the ubiquitous deco- 
herence effects disrupt the computation, and in general 
to define a more physical meaning to the complexity of 
quantum algorithms Q. The results published in time- 
optimal quantum control are numerous and it is certainly 
not the scope of the present paper to give an accurate or 
comprehensive review of this fascinating subject. Some 
literature may be found, e.g., in [To| . 

In a series of papers [ll[-[l4| we introduced a theo- 
retical framework for time-optimal quantum computing 
based on the action principle where the Hamiltonian is 
subject to a set of constraints (e.g., a finite energy, cer- 
tain qubit interactions are forbidden), and we called it the 
Quantum Brachistochrone (QB, [3I(). The QB was ap- 
plied to quantum state evolution of pure ll| and mixed 
states 12| , and to the optimal realization of unitary gates 
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. The predictions of the QB on the time complexity 
of generating unitaries with Hamiltonians which contains 
only one and two qubit interaction terms can be found 
in [Tj] . More recently, we discussed the problem of opti- 
mally generating entangling unitaries in a trilinear Ising 
chain subject to a coherent open loop control [lfj. As- 
suming a fixed energy available and that local operations 
on qubits are time-consuming, we found the laws of opti- 
mal control and the optimal minimal time for realizing a 
CNOT gate between indirectly coupled qubits. This time 
is shorter than previous results found in the literature so 
far [TEj . The importance of the study of quantum gates 
in physical Hamiltonians with indirect couplings is that 
non optimized procedures are usually time costly and eas- 
ily prone to decoherence effects and a degradation of the 
gate fidelity, both of which are critical for practical quan- 
tum computing. Interactions between indirectly coupled 
qubits via an intermediate qubit are also a typical sce- 
nario in a wide array of promising (scalable) experimental 
realizations of quantum information processing, e.g., in 
solid state architectures UM , crystal lattice architectures 
[r7| , superconducting [l8| and NMR technologies [l9[ . 

Here would like to elaborate upon [l(| on the prob- 
lem of time-optimal generating unitary gates in linear 
spin chains. The motivation is essentially twofold. On 
the one hand, we would like to introduce a general the- 
oretical framework for finding the time-optimal way of 
generating a unitary gate in the more realistic situation 
where the target can be reached within a finite, tolerable 
error, in other words with a fidelity close to but smaller 
than one. On the other hand, we would like to extend 
our previous analysis of Ising models to the case when 
the interaction couplings between spins are non e qual . 
Work related to our research can be found in (20j-[22l|. 
where analytical geometrical methods [23| and numerical 
algorithms [HI are used. We obtain new and interesting 
results for the minimal time duration of certain entan- 
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gling gates and propose some experiments to test our 
theoretical formulas. 

The paper is organized as follows. In Section II we 
briefly review the main features of the QBT formalism 
for the time-optimal synthesis of unitary quantum gates 
with a given final fidelity. In Section III we apply the 
formalism of the QBT and we introduce the main for- 
mulas for the problem of the efficient generation of the 
gate Ui 3 between the indirectly coupled boundary qubits 
1 and 3 of a three-linear Ising chain subject to a local, 
time consuming control on the intermediate qubit, when 
a finite energy is available. Section IV is devoted to the 
discussion of the case of a final perfect matching between 
the time-optimal unitary evolution and the target, while 
Section V analyzes the more general case when some er- 
rors are tolerated and the final fidelity is less than one. 
We show that in general the minimal time required to 
reach the target if smaller than previous estimates given 
in the literature. We also propose a series of possible 
experimental situations in which the results of the QBT 
formalism can be verified. Section VI is used to discuss 
the time-optimal generation with fixed fidelity of the gate 
CNOT(l, 3), which can be achieved using a slightly mod- 
ified interaction Hamiltonian plus the same local opera- 
tion on qubit 2 and the same energy available. Finally, 
Section VII is devoted to the summary and discussion of 
our results. The derivation of the main formulas exposed 
in the main Sections is outlined in Appendix A-D. 



II. THE QUANTUM BRACHISTOCHRONE 
WITH FIXED FIDELITY 

The goal is to minimize the time T it takes to perform 
a unitary operation under the condition that the fidelity 
/ between the time optimal unitary at time T and a cer- 
tain target is fixed. In other words, we want to realize a 
certain target unitary U /, but we tolerate that the target 
is reached only approximately, i.e. with a fidelity / < 1. 
We call this problem the Quantum Brachistochrone with 
Fidelity (QBF). We also assume that the Hamiltonian 
H is controllable and should obey some constraints, dic- 
tated either by theoretical conditions (e.g., only certain 
interactions among qubits are allowed) or by experimen- 
tal requirements (e.g., a finite energy), The unitary op- 
erator U(t) should then obey the Schrodinger equation. 
The QBF problem can be formulated in terms of the 
minimization of the following action: 

S(U, H, U(T);a 7 A, A, A,) := A[| (U(T), Uf) | 2 - (Nf) 2 } 

+ / ds [aN + L s + L c ] , 
Jo 

(1) 

L s :=(A,i^tf-aH), (2) 
/.«•: n^V/M//:. (3) 
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where (A,B) := Tr (A^B). The Hermitian operator A(t) 
and the real functions A and Xj(t) are Lagrange multi- 
pliers, while the Hamiltonian H{t) and the evolution op- 
erator U{t) are dynamical variables. The quantity \ogN 
represents the number of qubits. Also U (T) is a dynam- 
ical variable, since the goal in the QBF problem is to 
obtain the time-optimal trajectories reaching the target 
Uf only up to a fidelity /. The quantity a is the time 
cost. It plays the role of a "lapse" function [l2| connect- 
ing the physical time t := J a(s)ds to the parameter time 
s. 

The fidelity constraint is taken into account by the first 
term in the action S, and from variation with respect to 
A we obtain: 

|(^(T),^)| 2 = (V/) 2 . (4) 

Furthermore, variation of Ls by A gives the 
Schrodinger equation: 

i^=HU, or U(t)=Te~' l ti Hd \ (5) 
dt 

where T is the time ordered product. 

Variation of Lc by Xj leads to the constraints for H: 

/,(#) = 0. (6) 

In particular, the finite energy condition reads: 

UH):=\[Yt{H 2 )-Nu 2 ]=Q 1 (7) 

where a; is a constant. 

From the variation of S with respect to H we get: 

A = ^ + £^^ (8) 

where we have used ([7]). From the variation of S by 
a, upon using © and ([5]), we obtain the normalization 
condition: 

Tr(ZfA) = N, (9) 

while variation of S with respect to U(T) gives the fol- 
lowing final boundary conditions for A: 

A(T) = -2A lm[(U(T), U f ) U(T)u}}. (10) 

Finally, variation of 5" by U, use of eq. ^ and a simple 
algebra give the quantum brachistochrone equation: 

i^ = [H,A]. (11) 

The quantum brachistochrone together with the 
Schrodinger equation and the constraints define a two- 
boundary-value problem for the evolution of the unitary 
operator U(t) and the Lagrange multiplier A(t) with fixed 
initial U(t = 0) = 1 (where 1 is the identity matrix) and 
final conditions A(t = T) (cq. (|10p ). respectively. For a 
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given target gate Uf and fidelity /, the optimal Hamilto- 
nian H and the optimal time duration T can be found in 
the following way: a) solve the quantum brachistochronc 
(fTTj) backwards in time subject to the constraints © and 
the final boundary condition (fT0|) to obtain H opt (t); b) 
integrate the Schrodinger equation ([5]) forward in time 
with initial boundary condition U(Q) = 1 to get U op t(t); 
c) determine the integration constants in H opt (t) by im- 
posing the fidelity condition Q. 

III. THE QBF AND AN ISING HAMILTONIAN 

We now explicitly solve the QBF problem for the phys- 
ical system of three qubits (labeled by a superscript 
a € {1,2, 3}) arranged in a linear chain, interacting via an 
Ising Hamiltonian with interaction couplings J\ 2 , J23 and 
subject to a local and controllable magnetic field Bj(i) 
(i = x, y, z) acting on the intermediate qubit. More in 
details, we work with the Ising Hamiltonian: 

H{t) := J 12 a\a 2 z + J 23 a 2 z a 3 z + B{t) • a 2 , (12) 

where we have defined, e.g., o\g 2 '■= &i ® <Jj ® 1, the Oi 
are Pauli operators and N = 8. We simplify the notation 
introducing the ratio K := J23/ J 12 and rescaling the time 
as r := Ji 2 t, the energy as Cj := u>/ J12, the magnetic field 

as B(t) := B(t)/J 12 . 

Our goal is to time-optimally generate the symmetric 
cntangler gate U* 3 acting between the indirectly coupled 
qubits 1 and 3 (see eq. (4) in [IB]), i.e. the target is: 

Uf = U^:=e- i ^>" + ^ +a i\ (13) 

The form (|12[) of the physical Hamiltonian is guaran- 
teed by the operator (J5J: 

A(t) = \ h + E. tJ - fc \jk^]4 + E^W + dof] 

+ + ni<j\a) + Pn<y\(j% (14) 

where Ayfc(i), /3y (*) , Tfc (i) and are La- 

grange multipliers, and the indices {«,.?} € {2, y, z}. The 
finite energy condition ([7]) explicitly reads: 

B 2 = u% :=Qj 2 - (1 + A 2 ) = const, (15) 

and we must be sure that the available energy is enough 
so as to guarantee that uj 2 k > 0. 

The relevant set of equations associated to the quan- 
tum brachistochronc (|11[) with H and A respectively 
given by (|T2|) and (|14[) are presented as eqs. (|54l) - ([58|) 
in the Appendix A. In this Appendix we also show that 
Ao and v zz + Kp zz are integrals of the motion. There- 
fore, one can simplify the discussion by fixing the gauge 
freedom inherent to ((1) and choosing Ao = 1. From this 
and the quantum brachistochronc equation, it is then im- 
mediate to find out that B z = const and therefore, from 
the energy constraint (fT5"|) . that B 2 + B 2 := Bq = const. 



Exploiting these integrals of the motion and after some 
lengthy but elementary algebra, we find that the general 
and non trivial solution of the quantum brachistochrone 
equation (fTT]) is given by the Hamiltonian (|T^|) with the 
time-optimal magnetic field: 

(B cos 0(r)\ 
4 P t(r) = I B o sin0(r) j , (16) 

processing around the z-axis with the frequency f2, where 
9(t) := Or + 6q and Cl := O/J12 and 9q are integration 
constants. For later use, it turns convenient to use the 
energy constraint ([T5]) and to define (for 4> S [0,27r]): 

Bo := 0Jk cos 4>\ B z := lok sin (p. (17) 

The next step is to integrate the Schrodinger equation 
([S]) and to get the time-optimal evolution operator. The 
details of the calculation are shown in the Appendix B, 
and the result can be summarized as: 

f/ opt (T) = [a 13 (T)- l 6 13 (r).a 2 ], (18) 

where the operators a 13 (r) and b 13 (r) act in the Hilbert 
space of qubits 1 and 3 and are defined in eqs. ([77|) of 
the Appendix B. 

Finally, one has to impose the boundary conditions 
(ITU1) on A at the final time t* := J12T. Again, after some 
lengthy but simple algebra (for the explicit details the 
reader is referred to the Appendix C), it turns out that 
in the case of our 3-qubit model the relevant conditions 
can be summarized into the following equations: 

Or, = 2tott, (19) 

with m a non zero integer [32j and 

|M(t,)|=4/, (20) 
P(t.) = B(t*), (21) 

Q(t*) = §B(n), (22) 

where the functions M(t*), P(t*),Q(t*) and i?(r*) are 
defined in eqs. ([81~]) -([95" |) of the Appendix C. In other 
words, the QBF problem can be reduced to find the un- 
known integrals of the motion B z , Bq,Q,9q and t* from 
the constraints fl9"]) -([22 l) (via eqs. (JT3J) and ([67]l . see the 
Appendix C). 



IV. PERFECT MATCHING: / = 1 

We first consider the case of perfect matching between 
the time-optimal evolution operator t/ op t (■?"*) and the tar- 
get Uf — Ug 3 . In other words we assume that the fidelity 
is the maximum achievable, i.e. / = 1. In this case, it is 
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immediate to see from eqs. ([20 
the only possibility is to have: 



721-(|67P and ([STJ that 1 one gets: 



1,2,3,4, 



(23) 



where 712,^3,724 are positive integers such that n 2 = 
rii + 2q — 1, ri3 = 774 + 2s — 1 and 714 = n\ + 2r — 1, 
with q 1 r, s arbitrary integers. Then, from eqs. (|23[) and 
([72]) we find that all SiS are zero, and from eqs. (|86jl- 
([95]) we see that also P(t*) = Q(t*) = i?(r*) = 0, while 
M(r*) = 4(-l) ni . Therefore, eqs. (gni)-(E21) arc auto- 
matically satisfied, and one is left with the conditions 
(HU) and to be satisfied via eqs. ([73]) and flST]). This 
is essentially the same result which we had already found 
in our previous work [Io| for the case of equal couplings 
K = 1 (look at eqs. (43) and (63) in [Hj). In particular, 
following exactly the same methods of [10( , one first mul- 
tiplies eqs. ([73[) by r*, and then, using eqs. ([2"3"[) . inverts 
the former to find (with the help of ([P7j) ) the equivalent 
of eqs. (47)-(50) of 0, i.e.: 



T„ = 7T 




(So) 2 



8^ / fg^gZ 



K 

T 



2A/-IM!/, ^+.9 



.9+5- 



/- 



where we have introduced: 



/± := (n? + n|)±(nl + n^) 
.9± := (n? - "4) ± ( n l - "3) 



(24) 
(25) 

(26) 
(27) 



(28) 
(29) 



and their differences A/ := / + — /_, Ag := g + — g_. 
Moreover, the energy and the Ising interaction couplings 
are subject to the constraint: 



2K 



J K 



+ 4m m — 



9+9- 



V- 



- (1 + K% 



(30) 



the values of Ising couplings are also constrained by the 
relation K = g_ / g + [33j , and one has to guarantee that 
(-B0) 2 > and that uj 2 k > for the existence of the 
solutions [34j . 

Again, after a long but straightforward analysis (see 
[HJ), one can show that the perfect (/ = 1) time-optimal 
generation of the gate Z7-f 3 is possible for the values of 
/_ = 3 and ni even, and the solution coincides with 
that given in eqs. (51)-(54) of [l(J for the case of equal 
couplings. In particular, for m = 2 and 712 = 713 = 77,4 = 




and 



1 - 



— 777, + m 



(31) 
(32) 
(33) 
(34) 

(35) 



We stress that the duration time ([3"Tj) is shorter than the 
numerical result proposed in the earlier literature fl5| 
for the same trilinear Ising Hamiltonian and target as 
ours but for local controls which can selectively and in- 
stantaneously address the single qubits of the chain. In 
other words, as already cautioned by their authors, the 
quantum control procedure discussed in the paper [ljl is 
not time optimal (although it performs better than other 
standard procedures). 

For unequal couplings, i.e. K ^ 1, wc could not find 
solutions of the equations (IMl) - ([30[) which satisfy the con- 
ditions {Bo) 2 > 0, uj 2 k > and K = g-/g+ for /_ 
smaller than 7. The reason of this unexpected result for 
the K 1 case can be found by taking into considera- 
tion the more general ansatz of almost perfect matching 
between U op t{T*) and the target Uf, i.e. the ansatz of a 
fidelity / ^ 1. This is the topic of the next Section. 



V. ALMOST PERFECT MATCHING: / < 1 

We are now interested in finding the solution to the 
QBF problem in the case when the required final fidelity 
is not exactly one, but we allow for some error. In other 
words, we want to achieve a final fidelity: 



/ = 1 



(36) 



where |e| <C 1. By looking at the final boundary condi- 
tion (f2U[) and taking into account (|36[) , it is immediate to 
verify that the QBF problem can be solved by generaliz- 
ing the ansatz (|2"3"|) to: 



1,2,3,4, 



(37) 



where the real variables cr, 1, while the integers 71 j 
satisfy the same relations as those written below eq. (|23|) . 

Then, substituting ([37J into eqs. ([2"0" ]) -([2"2 ]l . wc obtain 
the following constraints for the o~iS: 



[bi 



b\o\ bnOi 



"4 



[63 



62CT2 &3C3 



"2 



(38) 



5 



(&x — sin <ft)— + (62 — sin0)— — (63 — sin 

m n 2 ' n 3 

,0"4 



, ^2 



^3 



condition A e x > or, in terms of the fidelity, using (|36|) 



and 



= (6 4 -sin0) — , (39) 
714 



(40) 



In principle, the procedure is now similar to that of the 
previous Section. In other words, one should invert eqs. 
(|73|) in order to express the unknown constants of the 
motion Bq,B z ,Q and t* in terms of the integers rii and 
the variables at, and then determine the former (con- 
stants) by minimizing r* over the latter (integers). How- 
ever, as it is shown in the Appendix D, only for values of 
the couplings ratio \K\ close to one it is possible to find 
time-optimal solutions. For pedagogical purposes, with- 
out losing generality we present the prototypical ansatz 
of positive K, i.e. K = 1 + 5K, with \8K\ < 1. The case 
of negative K can be dealt with along similar lines [35| . 

Then, performing a perturbative expansion in SK, af- 
ter some simple and tedious algebra (sec Appendix D) 
we obtain a self-consistent set of values for the variables 
<Ji (eqs. (|113[) - (|114|) of the Appendix D) which, for a 
given K close to one, leads to the following time-optimal 
constants of the unitary evolution: 



(41) 
(42) 
(43) 




(V2 + V3) 
+ (2V2 + V3) eK 



2(12^6 + 31) / A eK 



27 + 4^6 



where we have defined: 



<K 



SK — 



2(1 + 2^6) 



A, 



3tt 



27 + 4V6' 



(44) 

(45) 
(46) 



We point our that the optimal time duration (|4"Tj) is ap- 
proximately a function of A'" 1 / 2 , i.e. it decreases for in- 
creasing K. The optimal values of the integers are again 
n\ = 2, i%2 = 713 = n.4 = 1, for which /_ = 3, and in this 
case also the integer to can be optimized (see Appendix 
D) as to = 1. From the above equations, it is also clear 
that, in order for the time-optimal evolution to exist, i.e. 
for the values of r*, B z , Bq and fl to to be real, the ratio 
of the couplings in the Ising Hamiltonian must satisfy the 
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\SK\ < 3^VW- 



(47) 



This means that, the better we want to approximate the 
unitary target, i.e. the closer the fidelity should be to one, 
the smaller the deviation from one of the ratio between 
the Ising couplings that we can choose so that the pertur- 
bative expansion and the time-optimal evolution laws are 
self-consistent. Alternatively, if we have an experiment 
where K is given, we can read formula (|115[) as predicting 
that the distance of K from the unit value determines the 
maximal fidelity that we can hope to achieve (and impose 
in the QBF formalism). Furthermore, we have to choose 
an energy which satisfies the following constraint: 



(4V2 + 3^3) 




(48) 



We give a few examples of concrete situations where 
our formalism can operate and give the shortest time- 
duration of the quantum control in order to achieve the 
gate Ug 3 with fixed energy and fidelity. For instance, 
in the context of NMR |26|], we can think of an exper- 
iment where the Ising chain is the H — N — H part of 
the molecule of ethanamide [22[ , for which the Ising cou- 
plings arc J12 = J23 — 88.05 Hz^dt the F — F — F part 
of the trifluoroaniline molecule [23|, for which the Ising 
couplings are J12 = J23 — 20 Hz. For these molecules, 
K = 1 and the time optimal duration is given by eq. 
(f3Tj) of the previous Section. On the other hand, we can 
imagine an NMR experiment based on the H — H — H 
chain in the molecule of the l-chloro-2-nitro-benzene [28| . 
where the Ising couplings are J12 — 8 Hz and J23 ~ 7 
Hz, or an experiment based on the H — H — P chain in 
the molecule of 2-chloroethenylphoshonic acid [29| , with 
J12 ~ 9.1 Hz, J23 ~ 11.3 Hz. The former has a ratio 
K ~ 0.875, and the optimal fidelity reachable according 
to (|115|) is / max ~ 0.995, within a optimal time duration 
given by cq. (|4"Tj) . The latter has a ratio K ~ 0.805, 
and the optimal fidelity reachable according to (|115|) is 
/ max ~ 0.987. We notice that even taking the popular 
NMR model of the C — C — C chain the alanine molecule 
[30l] . with J12 ~ 54 Hz, J23 ~ 35 Hz and a relatively 
large SK ~ 0.35, the maximal achievable fidelity is still 
/ max ~ 0.957. It would be also interesting to apply the 
QBT formalism to experiments involving, e.g., supecon- 
ducting qubits. 



VI. CHANGING COORDINATES: CNOT(l,3) 

An analysis similar to that performed in the last Sec- 
tion can be done in the following two situations. On the 
one hand, the physical Hamiltonian (1121) may be available 
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in an experiment performed within the standard compu- 
tational basis (the Hilbert space spanned by the states 
{JO), |1)} <g> {|0), |1)} ® {|0), |1)}), but from the point of 
view of another lab the Hilbert space is seen as being 
spanned by the states {|0), |1)} <g> {|0), |1)} (8 {|+>, |-», 
where the basis for qubit 3 is {|+), |— )} (with |+) := 
W\0) = (|0) + |1))/V2; |-> := W\l) = (|0> - \l))/y/2 (W 
is the Walsh-Hadamard transform). Thus, the second lab 
sees an effective, rotated Hamiltonian given by: 

H'(t) = alal + Kalcxl + B(r) ■ a 2 . (49) 

Alternatively, one may simply assume that the Hamilto- 
nian (|49[) is available in an experiment performed within 
the standard computational basis for all the qubits. The 
Hamiltonians (jT2J) and (|4U)l are related by: 

H'{t) = V 3 H(t)V 3 , (50) 

where V 3 := l<g)l(g>W. Now, the goal is to time-optimally 
synthesize with fixed fidelity the gate CNOT(l,3): 

U f := CNOT(l,3) = e^iC^^"^-^), (51) 

As explained in our previous work [Toj . an analysis 
parallel to that of the previous Sections can be made. 
The relevant quantum brachistochronc equations are still 
given by (|54)) - (|58j) . with the only replacement pi Z — > pi X . 
Then the time-optimal evolution operator becomes: 

^opt(T) = V 3 U opt (r)V 3 , (52) 

where U op t(t) is given by (fT5)l . 

Also the new target CNOT(l,3) can be diagonalizcd 
(in the 1, 3 qubit subspace) as: 

U f = V 3 U'oV 3 = y 3 Diag(l, 1, 1, -l)^ 3 , (53) 

where the operator U' 1 ^ acts in the 1,3 qubit subspace. 
Then, following the same methods of the previous Sec- 
tion, it is easy to show that the only relevant changes are 
in a flip of the signs in front of the terms Si and S4 in 
eqs. (|8i~j) . (f86j) and (|95]l . Furthermore, in formula (|73|) 
we have choose the new positive integers n\ , i = 1, 2, 3, 4 
with n' 2 = ■n! 1 + 2q', n' 3 = n£ + 2s' and n' 4 = n[ + 2r' + 1, 
where q' , s' , r' are arbitrary integers. Considerations sim- 
ilar to those made in the previous Section finally lead to 
the expressions (I24ll - (l27|) . and to the minimization of the 
evolution time t* for /_ = 3 and n\ odd. In particular, 
for m — ri2 = U3 = 1 and = 2, the solution to the 
QBF problem with perfect matching (/ = 1) between 
U op t(T*) and the CNOT(l,3) target is given again, for 
equal couplings K = 1, by eqs. (f3"Tj) - (j35"j) (cf. eqs. (51)- 
(55) of [lOl )• For different couplings, the perturbative 
expansion in 5K leads instead to the time-optimal solu- 
tion (|41[) - (|4"8|) , except for a global change of sign for Cl 
and B z (see Appendix D). 



VII. DISCUSSION 

We have presented the general theoretical framework 
for determining the time-optimal way of generating a tar- 
get unitary quantum gate when the physical Hamiltonian 
is subject to certain constraints (at least a fixed energy 
available) and when some errors on the target are toler- 
ated, i.e. the fidelity is assumed to be fixed and smaller 
than one. We have then applied the formalism, which for 
brevity we called QBT, to the model study of a trilincar 
Ising chain with unequal Ising interaction couplings and 
where the intermediate qubit is subject to coherent con- 
trol via a local, but in principle time consuming, mag- 
netic field. For the case of a fixed fidelity equal to one 
(i.e. a perfect matching with the target), we explicitly 
gave analytical formulae of the time-optimal control laws 
and of the minimal time required to realize the entan- 
gling gate U} 3 between the indirectly coupled qubits of 
the chain when the Ising couplings are equal. The analyt- 
ical expression for the minimal time found is shorter than 
the numerical expression found (using standard geometri- 
cal quantum control methods for the same trilinear Ising 
Hamiltonian plus instantaneous individual control of the 
qubits) in [l5| . For the case of a fixed fidelity smaller than 
one, we performed an perturbative expansion around a 
ratio between the Ising couplings equal to one, and we 
gave analytical results for the time-optimal control laws 
and for the minimal time required to achieve J7 S 13 for un- 
equal couplings. Changing the computational basis for 
one of the qubits at the end of the chain, similar meth- 
ods can be used to analytically express the time-optimal 
control laws for generating the CNOT gate between the 
indirectly coupled qubits. All of the above results arc 
new, analytical and constitute an improvement with re- 
spect to the present literature for the values of the min- 
imal time required to generate certain unitary quantum 
gates between indirectly coupled qubits. We also pro- 
posed a series of experiments within the NMR paradigm 
to test our theoretical predictions. The main advantage 
inherent to the use of the QBF formalism appears to lay 
in the fact that (at least for these simple low dimensional 
models) it naturally allows for local quantum operations 
which in principle may have a non zero time cost, while 
the standard geometrical methods usually assume that 
local operations should be done instantaneously. In some 
sense, our ansatz is slightly more restricted than that of 
[l5| . since we also impose the finite energy condition eq. 
([7]). This, in turns, determines the appearance of the 
constraint (|35[) (or (|48p). which relates the allowed val- 
ues of the Ising couplings to the energy available in the 
experiment. As we already pointed out in (Toj . we con- 
jecture that this feature is the consequence of working 
with equality constraints on the Hamiltonian in the ac- 
tion principle, and we expect that it should be mitigated 
when such an assumption is relaxed and one considers in- 
equality constraints as well (for example, an energy avail- 
able bounded from above). We conjecture that also the 
difficulty in finding explicit results in the case of arbitrary 
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K far from one stems from the same technical point. We 
plan to extend our methods to the study of chains with 
more qubits and more complex topologies, to extend the 
allowed local controls to all qubits and to consider other 
quantum gates. Also in progress is some work related 
to the speed of quantum information transmission (e.g. 
transfer of coherence) along spin chains, and to the role 
of entanglement during the QBF evolution. 



VIII. APPENDIX 

A. QB and Integrals of Motion 

Comparing the coefficients of the generators of su(8) 
on both sides of (fTTj) , with A given by (fT4")) and H given by 
(I12p . we find that the relevant quantum brachistochronc 
equations are: 

(A SJ- = -2(u zy + Kp yz ), (54) 

(\ B y y = 2(v zx + Kp xz ), (55) 

(XoB z y = 0, (56) 
v zz = 2{B x v zy - B y v zx ), (57) 
p zz = 2(B x p yz - B y p xz ), (58) 

where Ao := J12A0 is a rescaled Lagrange multiplier. 
First of all, using the energy constraint equation ([7]) we 

note that (B) 2 = B x + B 2 + B 2 Z = const. Moreover, from 
dS]) , the quantum brachistochronc (fTTl) and (fT4"|) . we can 

compute = Tr[Ai2] = 8[w 2 A + {v zz + Kp zz )]. Then, 

taking the time derivative of (B) 2 , using eqs. (|54l) - ([55|) 
and the latter result, a simple algebra shows that Ao = d\ 
and is zz + K p zz = d 2 , where di,d 2 are constants. Finally, 
imposing the normalization condition one finds that 
ditb 2 + d 2 = 1. 



B. Time-optimal Evolution Operator 

We now proceed to integrate the Schrodinger eq. (O 
for U pt(T), given that H opt {T) is expressed via eqs. (O 
and (fT6|) . For this purpose, we exploit the following well 
known rotation formula for the Pauli matrices: 

e - i ° J i 1 ' y 'a x e iaJ i la * = cos6(r)a x + sm6{T)a yi (59) 



and we rewrite the time-optimal Hamiltonian as: 

■8(t) 2 -6(t) 2 



(60) 



where we have introduced the constant operator: 

H := B a 2 x + (a] + Kal + B z )a 2 z . (61) 
Furthermore, defining the transformed unitary operator: 

U( T ) := e^^Uir), (62) 



we easily check that, since U(t) should obey the 
Schrodinger equation ((5|), U{t) should also satisfy: 



if=HU 



with the time-independent Hamiltonian: 

O 



H := Hn- -at 



const. 



(63) 



(64) 



We note that the constant Hamiltonian H is diagonal in 
the 1,3 qubit subspace, i.e., 



const, 



(65) 



where we have introduced the operator B]^ acting in the 
1,3 qubit subspace: 



Brf := cj;fDiag[6i, b 2l 63,^4], 
which depends upon the constants: 



(66) 



sin < 



1 

0J K 



(5a-6u)(l + K) 



+ (6 i2 - 6 i3 )(l - K) 



(67) 



where i = 1,2,3,4 and Sij is the Kronecker symbol. 

Then, solving eq. (|63|) together with eq. (|65|) for U(t) 
and finally inverting ((62j) , it is easy to check that the 
time-optimal unitary operator U op t(T) evolves as: 



U opt (r) =e~ i * 



-iHr 



(68) 



In particular, the exponential of the constant Hamilto- 
nian appearing on the right hand side of eq. (|68|) is 
also diagonal (in the 1, 3 qubit subspace) and can be 
expanded as: 



e"^ = C 1 D 3 (r)- l S 1 D 3 (r)H, 
where we have introduced the diagonal operators: 



S 



Diag[si (r) , s 2 (r) , s 3 (t) , s 4 (r )] , 



Cif(r) := Diag[ci(r),c 2 (T),C3(r),c 4 (T)], 
which depend upon the functions of r: 
sin(o;iT) 



and the constants: 



Cj(r) := cos(wir), 



ujk\ cos 



bl 



(69) 



(70) 
(71) 



(72) 



(73) 
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Thus, inserting eq. (|69|) into eq. (|68|) . one obtains the 
more explicit expression (|18[) for the time-optimal evolu- 
tion of the unitary operator, where we have defined: 



where 



a (r) = cos 



Or 



C#(t)- sin 



Or 



6i d (r)=S ^(r)coB 



b?(T)=B Sg(T) sin 
Or 



Or 
~T 

Or 
~2~ 



bf (r) = sin 



C#(t) +cos 



Or 



B%Sh 6 (r), (74) 
(75) 
(76) 

Bi, 3 Si?(r). (77) 



The expression (fT5|) still depends upon the integration 
constants Bq,B z , O and #o an d the optimal duration time 
t * of the unitary evolution. These parameters are deter- 
mined by imposing the final fidelity constraint (j4]) and 
the final boundary condition (fTU| . 



C. Final Boundary Conditions on A 

The boundary conditions for A at the final time r* 
J12T in our 3-qubit model explicitly read: 

A(t*) - 2A C/^ 3 6 13 (t„) ■ a 2 Tr[C/ opt (r,)C/)]. 



(78) 



Let us now compute these conditions. First, we note that 
our final target (|13p is diagonal in the 1, 3 qubit subspace 
and can be written as: 

U f = e i *Uh 3 , :=Diag(-l, 1,1,1). (79) 

Then, using (JTSJ) and ([79")l . we can write: 



Tv[U opt (T*)U f ] = 2 (cos 



Or* 


M — sin 


Or* 


2 




2 



(80) 

where we have defined: 

M(t.) := TrtCi, 3 ^ 3 ] = Cl + c 2 + c 3 - c 4 , (81) 

JV(t„) := Tr^ 3 ^ 3 ^ 3 ] = Mi + &a*a + Ms - M 4 . 

(82) 

Now, inserting (|50|) into ([75)) , from the a 2 , er 2 , , er 2 terms 
of A(r*), we obtain the final boundary conditions: 



cos 0(t„) 



AP(t.) 
2wx 



cos 



AP(r») . 

sin 

2wk 



Or* 
Or* 



Tr[C7 opt (r*)C/J], 

(83) 

Tr [[/ opt (t* )[/}], (84) 



-I cos 



Or* 


N + sin 


Or* 




2 




2 



AI 



Tr[U opt (T*)U f ], 



(85) 



P(t») := Tt[S%Ug] = 81 + s 2 + s 3 - fl4 . 



(86) 

Upon using (JHOJ) into ([Ml) and (J5JJ we find that these 
can be nontrivially [36( | satisfied only if: 



(87) 



tan 0(t») = tan 

whose solution is given by: 

Or, = 2to7t, 



where m is a non zero integer [371 ] . Using the latter for- 
mula we can compactly rewrite equations (|83[) -(|85j) as: 



AM(t,)P(t ( ) = wjr, (89) 
AA/(t*)^(t») = (90) 
Moreover, from the fidelity eq. ((4]) and the definition: 

Tr[U opt (T*)U}} := | Tr[[/ opt (T»)[/t]| e - 4 * = 8/ e" 4X , 

(91) 

(92) 
and eq. ([20| . 

Let us now impose the normalization condition ^ at 
the final time r». Multiplying (fTOj) by H(t*) and using 
([9|); we can express the Lagrange multiplier A as: 



and from (|80|) and (|9Tj) . we obtain: 

e 8 * = (-l) m sign[M(r*)] 



A = -{2 ^(e^TrlC/optCT^C/jHopt^)])}- 1 . 

After some lengthy but simple algebra, from (fT3 
(fI5]l and dHFJ one can compute: 



(93) 



Tr[U opt {T*)u}H opt {T*)} = 2i{-l) m+ W K R{T*), (94) 



where we have defined: 



R{i 



cos 

2 



2 sin < 



cos 



cos 



cos" 



+ — (fo - h + 2sin< 
- ^(i>2 - h - 2sin< 
— —(61 — &4 — 2 sin< 



}S 4 . 



(95) 



Finally, substituting ((92|) and ((94]) into ([93]). we get: 
sign[M(r*)] 



A 



4fu 2 K R(T*) 



(96) 



Thus, if we can determine the values of M(r*) and i?(r»), 
from the latter equation we have determined the optimal 
value of the Lagrange multiplier A for a given fidelity 
Furthermore, using (|96[) the multiplier A can be elimi- 
nated from ([M)) - (f9"0"]) which, after some simple manipula- 
tions, become eqs. (|2"Tj) - ([2"2"j) . where 

Q(t») := u K [(l + K )( s i + s i) + (1 - K )( s 2 - ss)]- 

(97) 
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D. Perturbative Expansion around K = 1 

Wc first compute the function /_ for the case when the 
target is the entangler gate U} 3 , i.e. when the integers are 
n\,ri2 = ni+2g— l,n3 = ni+2s — 1 and 714 = m+2r — 1. 
From its definition, eq. (|2"9")l , we get: 



/_ = 2m[2(r-flf-a) + l] 

+ 4[r(r - 1) - q(q - 1) - s(s - 1)] 



(98) 



Now let us assume that the Ising couplings are not equal, 
i.e. that \K\ 7^ 1, and let us compute (wit*) 2 , (W2T*) 2 , 
(W3T*) 2 and (aj4"rv) 2 . From the definitions (|57|) . ([75)1 and 
cqs. (fTTj|) and (j3"T)) we obtain: 



(cj^t*) 2 + 2(lj kt*) sin</>[(l + K)t* — mn} 

+ [(1 + K)t* - mn} 2 = (win + ai) 2 , 
(ujkt*) 2 + 2(ujkt*) sin</>[(l — K)t* — mn] 

+ [(1 - K)t* - mn} 2 = (nn 2 + cr 2 ) 2 , 
(uj k t*) 2 — 2(<jj k t*) sin</>[(l — K)t* + mn] 

+ [(1 - K)t* + mn} 2 = (nn 3 + a 3 ) 2 , 
(ujkt*) 2 — 2{ujkt*) sin</>[(l + K)t* + mn} 

+ [(1 + K)t* + m7r] 2 = (7rn 4 + cr 4 ) 2 - 



(99) 
(100) 
(101) 
(102) 



Then, summing to f|102[) and subtracting from this 
(jlOOP and (jlOip , we obtain the duration time of the uni- 
tary evolution: 



2/f 



(103) 



where we have introduced (for notational consistency 
with (|2"9")l ) /J := {ri\(j\ + 71404) ± (n 2 03 + n 3 cr 3 ), and 
we discarded contributions of O(o 2 ). Subtracting (jlOip 
from pOOp . we obtain the constant component B z of the 
magnetic field: 



u>k sin ( 



n 2 Ag 



[1 - K)t. 



2Ag c 



nAg 



(104) 



where g a ± := (nioi - n 4 cr 4 ) ± (n 2 cr 3 - n 3 cr 3 ), Ag a := 
g'L — g°_ (contributions of 0(a 2 ) discarded). From the 
sum of dS9j), (fT00|) . (fTDTjl and (fT02|) . wc get instead a 
constraint on Cj 2 (and therefore, via p~5|) . on w|-): 



7T 

ir 2 



/+ + 4m 2 



mnAg + 2fl 
2t~ n 



(105) 



Finally, from the sum of (|99|) and (|10ip . to which we 
subtract (jlOOp and (|102p . we obtain the constraint: 



(l-K)(g + +g_)-(l + K)Ag 

= -[{l + K)Ag-{l-K){g^ 
n n 



9-)}- (106) 



(modulo 0(a 2 ) contributions). The last equation ex- 
plains the origin of the unexpected result of Section IV. 
In fact, for \K\ ^ 1, the only way for which cq. (j 106[) can 
hold to zero order in a is that the constraint K = </_ / g + 
must hold. When one is looking for the minimum (odd 
integer) value of of the function /_ (i.e., from f)103[) the 
minimum time duration r* of the unitary evolution) this 
constraint together with those coming from (Bo) 2 > 
and uj 2 k > (see the discussion below cqs. (f24"]) - (j30|) ) 
must be simultaneously satisfied by the integer numbers 
m, q,r and s. This is a very strong set of constraints. 
No surprise that, at least for /_ up to 7 (we did not try 
higher, unrealistic values of /_ and, therefore, of r»), only 
the choice of n 2 = n 3 is allowed, thus enforcing g+ = gr- 
and K = 1. The only possible way out of this "empasse" 
is to assume that we can work with \K\ close to one and 
hope that the right hand side and the left hand side of 
eq. (|106p are of the same order. 

Let us then assume, for example, that K = 1 + 6K, 
where we take \SK\ <C 1. First of all, performing an 
expansion m SK of the condition fl06]) . we indeed find 
that this separates into two conditions, coming from the 
zero order and the first order terms in SK and cr,, i.e., 
respectively: 



m = m, 



712(0-2 - 0-3) 



n r 2 



n\)8K. 



(107) 
(108) 



That is to say, SK must be of the order of o^, i.e., from 
(UOJ), we must have SK iS e. 

Then, we perform a perturbative expansion in SK of 
the conditions on the variables Oi, i.e. an expansion of 
eqs. (|38l) - ((39|) . Exploiting (|67| we obtain, up to the 
lowest order in both SK and o;\ 



u>k sin ( 



"2 



1L 



u>K sm < 



(- 


a 4 \ 


n / 


v m 


n 4 J 


~ 4 [ 



"1 



(72 

ri2 



' 2 

m 

n 3 



- 2 



1± 

»4 



"4 



(109) 



(110) 



We can now solve eqs. (|109l) - (|110p together with ([Hi]). 
Using the lowest order approximations to the time- 
optimal evolution (eqs.([31" ]) - ([33"]) . or cqs. (fT03 l) - (fT04| . 
with n\ — 2,ri 2 = 71,3 = n& — 1 and /_ =3), we can 
find that w^-sm</>— Cl/2 ~ — 1. Therefore, eqs. (| 1 09[> - 
(|TT0| give: 



0"2 = - 



^ ) (j 1 - *L6K 
ml 4 



03 



V6 N 



cri 



3tt 



5K 



CT4 = 



(Ill) 

(112) 
(113) 
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Finally, substituting (|113[) into (|40|) . wc obtain: 



ai=± 



\ 



2A, 



K 



l+yjf) ,13 



(114) 



where A e x := e 2 — ^-(Jif) 2 . Here we notice that, in 
order for the time-optimal evolution to exist, i.e. for the 
Cj to be real, we have a more precise condition that must 
be satisfied by the couplings in the Ising Hamiltonian, 
i.e., we must have that A e x > or, in terms of the 
fidelity, using eq. (|36|) : 



1G 



\aki < -yw. 



(115) 



We can now proceed in determining the perturbative ex- 
pansions in SK of the time-optimal solutions (|103p - (|106[) 
for /_ =3, which explicitly read (to the lowest order in 
SK and <x,-): 




(116) 
(117) 
(118) 



U ~ m\ — 1 H 

V 3 V 2 



where wc have defined: 



S Ka := SK - 



2fZ 
3tt 



(119) 



Finally, the value of the integer m is also fixed by mini- 
mizing t* in (|116p . i.e. by minimizing 



1 



2\/6 N 



0"! 



(120) 



with <j\ given by f|114[) . with respect to m (and the sign 
of a \ in (| 1 14[) ^ . The result is m = 1 and sign(<7i) = +, 
for which f°_ = -(1 + 2V6)[A eK /(27 + 4 V / 6)] 1/2 . This 
gives the final formulas ([4l"]) - (|43]) of Section V. 

A similar computation can be done for the case of 
the target CNOT(l,3), cq. flST), and for the "rotated" 
Hamiltonian (|49l) . discussed in Section VI. In particular, 
i' 2 = n[ + 2q',n' 3 = n[ + 2s',n' 4 = 
n^n^n's and n 4 positive integers, 
and q',s' and r' integers. One obtains, for m = 2 and 
112 = «3 = «4 = 1, the same time-optimal solutions (|41jl - 
(|4"8")) . except for a change everywhere in in the sign in 
front of odd powers of the integer m, due to the fact that 
now m = — 1 is the value minimizing t*. 



one has to choose i 
n\ + 2r' + 1, with 



ACKNOWLEDGEMENTS 

This research was partially supported by the MEXT 
of Japan, under grant No. (T.K.) and by the project 
Rientro dci Ccrvclli of the MIUR of Italy (A.C.). 



[1] M. Shapiro and P. Brumer, Principles of the Quan- 
tum Control of Molecular Processes (Wiley, New Jersey, 
2003). 

[2] M.A. Nielsen and I.L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

[3] C. Brif., R. Chakrabarti and H. Rabitz, New J. Phys. 12, 
075008 (2010). 

[4] J. Cavanagh, W.J. Fairbrother, A.G. Palmer and N.J. 
Skelton, Protein NMR Spectroscopy: Principles and 
Practice (Academic Press, San Diego, 1996). 

[5] N. Margolus andL.B. Levitin, PhysicaY>120, 188 (1998). 

[6] V. Giovannetti, S. Lloyd and L. Maccone, Europhys. Lett. 
62, 615 (2003); ibidem, Phys. Rev. A67, 052109 (2003). 

[7] M. Murphy, S. Montangero, V. Giovannetti and T. 
Calarco, Phys. Rev. A82, 022318 (2010). 

[8] M.A. Nielsen, M. Dowling, M. Gu and A. Doherty, Sci- 
ence 311, 1133 (2006); id., Phys. Rev. A73, 062323 
(2006). 

[9] T. Schulte-Herbriiggen, A. Sporl, N. Khaneja and S.J. 
Glaser, Phys. Rev. A72, 042331 (2005). 
[10] A. Carlini, A. Hosoya, T. Koike and Y. Okudaira, J. 
Phys. A44, 145302 (2011). 



[11] A. Carlini, A. Hosoya, T. Koike and Y. Okudaira, Phys. 

Rev. Lett. 96, 060503 (2006). 
[12] A. Carlini, A. Hosoya, T. Koike and Y. Okudaira, J. 

Phys. A41, 045303 (2008). 
[13] A. Carlini, A. Hosoya, T. Koike and Y. Okudaira, Phys. 

Rev. A75, 042308 (2007). 
[14] T. Koike and Y. Okudaira, Phys. Rev. A82, 042305 

(2010). 

[15] N. Khaneja, B. Heitmann, A. Sporl, H. Yuan, T. Schulte- 
Herbriiggen, and S.J. Glaser, Phys. Rev. A75, 012322 
(2007). 

[16] B.E. Kane, Nature 393, 133 (1998). 

[17] F. Yamaguchi, Y. Yamamoto, Appl. Phys. A: Mater. Set. 
Process. 68, 1 (1999). 

[18] J. Clarke and F. Wilhelm, Nature 453, 1031 (2008). 

[19] R.R. Ernst, G. Bodenhausen and A. Wokaun, Principles 
of Nuclear Magnetic Resonance in One and Two Dimen- 
sions (Clarendon, Oxford, 1987). 

[20] H. Yuan, R. Zeier and N. Khaneja, Phys. Rev. A77, 
032340 (2008). 

[21] H. Yuan, S.J. Glaser and N. Khaneja, Phys. Rev. A76, 
012316 (2007). 

[22] M. Nimbalkar, R. Zeier, J.L. Neves, S.B. Elavarasi, H. 



11 



Yuan, N. Khaneja, K. Dorai and S.J. Glaser, Phys. Rev. 

A85, 012325 (2012). 
[23] N. Khaneja and S.J. Glaser, Chem. Phys. 267, 11 (2001); 

N. Khaneja, R. Brockett and S.J. Glaser, Phys. Rev. 

A63, 032308 (2001). 
[24] N. Khaneja, T.O. Reiss, C. Kehlet, T. Schulte- 

Herbriiggen and S.J. Glaser J. Mag. Res. 172, 296 

(2005). 

[25] B. Bonnard, S.J. Glaser and D. Sugny, Adv. Math. Phys. 

2012, 857493 (2012). 
[26] J. A. Jones, arXiv:quant-ph/1011.1382. 
[27] O. Mangold, A. Heidebrecht and M. Mehring, Phys. Rev. 

A70, 042307 (2004). 
[28] D.G. Cory, M.D. Price and T.F. Havel, Physica D120, 

82 (1998). 

[29] H.K. Cummins, C. Jones, A. Furze, N.F. Soffe, M. Mosca, 
J.M. Peach and J. A. Jones, Phys. Rev. Lett. 88, 187901 
(2002). 



[30] D.G. Cory, M.D. Price, W. Maas, E. Knill, R. Laflamme, 
W.H. Zurek, T.F. Havel, and S.S. Somaroo, Phys. Rev. 
Lett. 81, 2152 (1998). 

[31] From the Greek "^pa^Kjro^", i.e., fast, and "\povoQ\ 
i.e time. 

[32] If m — 0, either T = or O = 0, which are both trivial 
cases. 

[33] In fact, when the couplings are equal, then g+ = g~ 
and coincide with /o of eq. (46) in [10]. Then from 112911 
ni — 713 so that also the f± coincide with eq. (45) of [10]. 

[34] This, in particular, requires that sign(/_) = sign(if) = 
sign(#+3-)- 

[35] In this case, for the reality of the duration time, eq. (|31l) . 

one has to work with /_ negative. 
[36] The trivial solutions are Bo or B z equal to zero. 
[37] If m — 0, either T = or O = 0, which are both trivial 

cases. 



